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Abstract. We investigate the behavior of the Green functions of Schrodinger 
operators near the diagonal. The only non-trivial cases, where the on-diagonal 
psj ■ singularities are non-zero and do not depend on the spectral parameter, are two and 

three dimensions. In the case of two dimensions, we show that the singularity is 
independent of both the scalar and the gauge potentials. In dimension three, we 



00 ' obtain conditions for preserving the singularity under perturbations by non-regular 

. potentials. Some examples illustrating dependence of the singularity on general 

I I scalar and gauge potentials are presented. 

o 

i!^ ' 1 Introduction 

f~| ■ Singularities of the Green functions of the quantum-mechanical operators play a crucial role in 

"j^ I many branches of theoretical and mathematical physics, from which one should mention first the 

renormalization procedure of the quantum field theory [1-3]. From the point of view of the high- 
derivative quantum gravity, the corresponding problem was considered e.g. in [4]. In particular, 
in the case of non-minimal coupling of quantum matter to the gravitational background with 
^ ' conical singularities, an operator of the form H = —A + U arises on a Riemannian manifold X. 

■ Here A is the Laplace-Beltrami operator on X and U represents the non-minimal coupling term 

S^TZ with the Ricci scalar TZ. The scalar curvature possesses a distributional behavior at conical 
singularities [5], TZ = T^rog + 47r(l — a)5M , where 5m is a Dirac (5-like potential supported by a 
sub-manifold M G X and 27r(l — a) is the angle deficit. As a result, an operator 

Hm = -A + V + a5M (1.1) 

arises with the coupling constant a = 47r(l— a)^ characterizing the interaction with a background 
field concentrated on M. Operators of such form appear in the investigation of scalar fields with 
non-minimal coupling on the cosmic string background, in the Euclidean approach to the black 
hole thermodynamics, in the study of the particle scattering at the Planck scale (see [5] and 
references therein). Moreover, in the context of the scattering theory, the potential V can have 
singularity (e.g. of the Coulomb type) even in the case of a flat manifold X. 

We are interesting here in the singular term 6m concentrated on a zero-dimensional sub- 
manifold M; this case covers not only quantum fields with point interactions, but also the 
case when M has a cartesian complement in X: X = Y x M. If M is a uniformly dis- 
crete subset of X and dimX < 3, then the Green function GMix,y',C) of H can be obtained 
through the Krein resolvent formula in terms of the Green function G{x, y; for the operator 
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H = —A + V [6]. An important ingredient of this formula is the so-called "Krein Q-matrix" 
(a kind of the Dirichlet-to-Neumann map) Qmn{z), m,n ^ M. To define the diagonal elements 
of Q for dimX > 1, a renormalization procedure is needed. For smooth V, the renormalized 
Green function G'^°^{x,yX), which must be continuous in the whole X x X, is defined as 

G''-(x,y;C)=G(x,y;C)-S(x,y), (1.2) 

where the "standard singularity" S has the form S{x,y) = — — log d{x , y) if dimX = 2, and 

S(x,y) = - — r- if dimX = 3 (here d(x,y) is the geodesic distance on X). Now one can put 

47ra(x, y) 

QmmiC) = G^^'^{m,m; C). The corresponding renormalization procedure in the Euclidean case 
is known long ago, see e.g. [7] and [8] for the history and the quantum mechanical treatment. 
It is important to note that usually one obtains S{x, y) by a momentum cutoff (an ultraviolet 
regularization procedure); the result is equivalent to that obtained with the help of a dimensional 
regularization. In the case of brane coupling to gravity or to a gauge field it is necessary to use 
a dimensional regularization [9]. It is worthy to add that the strict mathematical treatment of 
the operators (|1.1|) has its origins in the article [10] by F. Berezin and L. Faddeev. 

In the case dimX > 4 there is no regularization procedure involving a singularity indepen- 
dent of the energy parameter ^ (see Example [7] below) . Moreover, if V has a Coulomb-like 

singularity or if an interaction with a gauge field is present, then the function 5 in p.2|) is 

1 

different from the standard one, i.e., S{x,y) ^ (see Examples 1121 and 1141 below) : sim- 

ATrd{x, y) 

ilar phenomena related to propagation of waves in strongly inhomogeneous media have been 
studied recently in [11]. The main goal of our paper is to investigate the situation in detail. 
We show that in dimension two the singularity of G has the standard form even in the presence 
of an additional U(l)-gauge potential (Theorem I15() . On the other hand, in dimension three, 
S depends on V modulo a Lebesgue class of functions on X (see Theorem ITH]) and is defined 
up only to a continuous additive term (the situation here is completely similar to that for the 
Krein Q-functions: they are defined up to an additive constant). The concrete value of this term 
is subject of analysis of a given physical problem and is out of the scope of the present work. 
We mention only that a possible way to fix the corresponding additive constant is to compare 
the integrated density of states with the trace of C^*^". It is worthy to note that the Green 
function for operators of the form Hl.l|) on a three-dimensional Riemannian manifold was used 
recently for simulating the confinement potential of a quantum dot in [12]. Earlier the defects in 
solids are investigated by methods of quantum gravity in [13]. New technologies of manufactur- 
ing two-dimensional nanostructures with non-trivial geometry [14, 15] caused the appearance of 
mathematical models of such structures where, in particular, the Hamiltonian has the form (|1.1() 
with the 5-term simulating the potential of a short range impurity [16]. If the nanostructure is 
displaced in a magnetic field we must replace A in (jl.lj) by the Bochner Laplacian. In this case 
the properties of the Green function G are needed for investigation of explicitly solvable models 
of the geometric scattering theory [17] or spectral theory of periodic hybrid manifolds [18]. 

At last but not at least we stress that our main results are new even for the case of Euclidean 
spaces X = W^. 



2 Definitions and preliminaries 

Throughout the paper we denote by X a complete connected Riemannian manifold of bounded 
geometry, which means that the injectivity radius rinj of X is positive and every covariant 
derivative of the Riemann curvature tensor is bounded. Examples are provided by homogeneous 
spaces with invariant metrics (in particular, Euclidean spaces), compact Riemannian manifolds 
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and their covering manifolds; for discussion of various properties of such manifolds in the context 
of differential operators we refer to [19]. The dimension of X we denote by u; the geodesic 
distance between x,y £ X will be denoted by d{x, y). For x £ X and r > we use the notation 
B{x,r) = {y £ X : d{x,y) < r}; through the paper, we suppose r < rinj for radiuses r of all 
considered sufficiently small balls. For a measurable function / on X, we denote by ||/||p the 
L^-norm of /. If is a bounded operator from U'{X) to L'^{X), 1 < p,q < oo, then its norm 
will be denoted by 

Let A = Aj dx^ be a 1-form on X, for simplicity we suppose here Aj £ C°°{X). The functions 
Aj can be considered as the components of the vector potential of a magnetic field on X. On the 
other hand, A defines a connection in the trivial line bundle X x C ^ X , = du + iuA; 
by —Aa = V^Va we denote the corresponding Bochner Laplacian. In addition, we consider a 
real- valued scalar potential U of an electric field on X. This potential will be assumed to satisfy 
the following conditions: 

n 

U+ := max(;7,0) £ L^^^iX), U. := max(-C/,0) £^LP^{X), 

1=1 

2 < Pi < oo if I' < 3, Z//2 < Pi < cxD if 1/ > 4, < i < n; 

we stress that pi as well as n are not fixed and depend on U. The class of such potentials will be 
denoted by V{X). Below we will need an approximation of singular potentials by smooth ones; 
for this purpose the following lemma is useful. 

Lemma 1. Let f £ L^^^{X), where 1 < p < oo, and / > 0. Then there is g £ C^{X) such that 
g>0 and f -g £ L''{X) for alll<q<p. 

Proof. Fix a £ X and for integers n, n > 1, denote Yn = B{a,n) \ B{a,n — 1). Fix a real 
sequence a^, > such that X] — ^ denote by fn the restriction of / to the set Yn- 
Since the measure of Yn is finite, for every n we can find a function gn, gn S C^{X), such 
that gn > 0, supp(5(„) C and max(||/„ - gn\\p, \\fn - ffnlli) < a„. Since the family (y„) is 
locally finite, the point-wise sum g = Yldn exists and g £ C^{X). It is clear that g > and 
max(||/ - g\\p , ||/ - 5II1) < 1, i.e., f-g£ LP{X) nL^X). □ □ 

We denote by Ha,u the operator acting on functions (p £ C^{X) by the rule Ha,u4' = 
—AA(j) + U(p. This operator is essentially self-adjoint in L^(X) and semibounded below [20]; its 
closure will be also denoted by Ha,u- By spec{HA,u) we denote the spectrum of Ha,u and by 
res{HA,u) the set of regular points: Tes{HA,u) = C \ spec{HA,u)- Let us denote the resolvent of 
Ha,u by Ra,u{C), i.e. i?A,(/(C) = {Ha,u - C)"'- 

Here we introduce two classes of integral kernels used in the paper. First class, /Ccont(p)) 
1 < p < 00, consists of all continuous on X x X functions K{x,y) satisfying for any r > the 
condition 

[K\p^r ■■= max sup ess^.^^ \\xx\B{x,r) K{x, ■)\\^ ,sup esSy^^x \\xx\B(y,r) ^(-^ J')||p) < 00 ' 

(2.1) 

where XA stands for the characteristic function of the set A C X. The second class, IC {a,p), < 
a < z/, 1 < p < 00, consists of all measurable functions K on X x X obeying the condition (|2.1j) 
and 

|-f^(2;,y)| < cmax(l, for a constant c = c{K) > 0. (2.2) 

We put K, cont{oi,p) := IC {a,p)nC{X xX\D), where D is the diagonal {{x,y) £ XxX : x = y}. 

The above introduced classes of integral kernels are important due to their relations to the 
properties of the resolvents RA,uiC)] these relationships are stated in the following theorem 
which is our starting point (see [20] for the proof). 
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Theorem 2. For any C, G ies{HA,u) the resolvent Ra,u{C) has an integral kernel GA,u{x,yX), 
the Green function, which belongs to /C cont(A, 9), where q, i <q < 00, is arbitrary, and X = v — 2 
for > 2, A G (0, is arbitrary for = 2, A = for v = 1; moreover, Ga,u is continuous in 
XxXforu = l. 

We should point out that the Green function of a Schrodinger operator can violate the 
conditions ()2.1|) and ()2.2() . if the potential U is not from the class V{X). Even the decay of the 
Green function for large distances between x and y (the off-diagonal behavior) can be different 
from the "standard" exponential one coming from the comparison with the Laplacian; a good 
example is delivered by the one-dimensional inverse harmonic oscillator, whose Green function 
has only a polynomial decay at infinity (see Appendix^. 

Our further calculations will involve a couple of operations with integral kernels introduced 
above; here we collect some useful estimates which will be used very intensively. 

The well-known Gelfand-Dunford-Pettis theorem claims that if is a bounded operator 
from U'{X) to L°°{X) with some p, 1 < p < 00, then it is an integral operator and its kernel 
K{x, y) satisfies the estimate 

sup eSS^gj^ 11^(2;, Ollg < 00, gr = (1 -p-l)-l. (2.3) 

Conversely, if a kernel K{x, y) satisfies 1)2. 3() . then it is an integral kernel of a bounded operator 
from LP{X) to L~(X). 

Lemma 3. Let Kj : L'^j(X) L°^'{X), 1 < qj < oo, be bounded linear operators with integral 
kernels Kj{x,y), j = 1,2, and W G L^'^{X), then for a.e {x,y) G X x X the integral J{x,y) = 

Ki{x, z)W{z)K2{z,y)dz exists and J{x,y) is an integral kernel of the operator K1WK2. 

X 

Proof. The operator K1WK2 is bounded from L'^'^{X) to L°°{X), therefore, it is an integral 
operator. Let / G L92(X) n C{X) such that f{x) > for ah x G X. Then there holds 



K^WK2f{x) = / K^{x,z)W{z) / K2iz,y)f{y)dydz. (2.4) 
Jx Jx 

From the other side, according to the estimates 1)2. 3|) for Ki and K2, there holds 

/ \K2{;y)f{y)\dy€L^{X), W{-) [ K2i;y)fiy)dy€L'^^{X), 
Jx Jx 

hence, 

^ \Ki{x,z)\ (^\W{z)\ \K2{z,y)f{y)\dy^dz < 00 . 

By the Fubini 

^ (^j^\Ki{x,z)W{z)K2{z,y)\dz^ f{y)dy < 00 , 

and since f{x) > 0, the inner integral exists for a.e. (x,y) G X x X. 

Let now / be an arbitrary function from L'^'^(X). Repeating the arguments above, we get 

KiWK2f{x) = (^j^Ki{x,z)W{z)K2{z,y)dz) f{y)dy (2.5) 

for a.e x £ X. Therefore J is an integral kernel for KiWK2- □ 
We will often use the estimate given by the lemma below (cf. [20]): 
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Lemma 4. There exists > such that for any a, r with < r < tq, < a < i/, and a,x ^ X 
there holds 

I < cr'^"" (2.6) 

with some c > depending only on a. 

Our next auxiliary result is the following lemma. 

Lemma 5. Let K £ IC {a,p), I < p < oo, pa < v, and 1/p + 1/q = 1, then K is an integral 
kernel of a hounded operator from L'^{X) to L°°{X). 

Proof. According to the Gelfand-Dunford-Pettis theorem we must prove 

sup ess^^x \K {x, y)\^ dy < oo. 
Jx 

Fix r, < r < ro, and for x £ X expand the integral into two parts: 

K{x,y)\''dy= [ \K{x,y)\''dy+ [ \K{x,y)\''dy . 

X ' JB(x,r) Jx\B(x,r) 

The first term is estimated by Lemma 0] and the second one is maj orated by [i^Jp^r- Q 

Lemma 6. Let three measurable functions Ki(x,y), K2{x,y) and W{x) be given, where x,y £ 

X. Denote F{x,y, z) := Ki{x, z)W{z)K2{z,y), and if the integral / F{x,y, z) dz exists, denote 

Jx 

it by J{x,y). 

(A) Let Kj G fC contictj,Pj), j = 1,2, and W £ LP{X), such that l/pi + l/p2 + 1/p = 1 o-nd 
p > v/ {y — max(ai, 02)) . Then F{x, y, •) £ L^{X) for x ^ y, hence J is well defined. Moreover, 
J G ^cont(a, 00); where a = max {p'{ai + 02) — ^^,0) with 1/p + 1/p' = 1, if p'{ai + 02) 7^ i^, 
and a is an arbitrary number from (0, v) otherwise. 

(B) Let the conditions of the item (A) be satisfied. Assume additionally that ai + 02 < f 
and W £ L1^^{X) with q> v/{v — a\ — 02)- Then F{x,y, ■) £ L^{X) for any x,y £ X and 
J £ C{X X X). 

(C) Let W £ LP{X), and Ki £ /Ccont(pi), e K^contia,P2) or Ki £ /C cont(a,Pi), e 
/Ccont(P2)- Assume additionally that l/p + l/pi + l/p2 = 1 and p > u/{u — a). Then F{x,y, •) £ 
L^{X) for any x,y £ X, and J £ C{X x X). 

Proof. The proof of the items (A) and (B) is given in [20] . 

(C) We give a proof for the case Ki £ /Ccont(pi) and K2 £ cont(cK)P2); the second case can 
be considered exactly in the same way. 

Let x,y £ X; we show first that F{x, y, •) £ L^{X). Let r > 0, then for z £ B{y, r) we have 

\F{x,y,z)\ < cki{x,y)W{z)d{y,z)^°' , ki{x,y) := sup Ki{x,z)<Qo, c > 0, (2.7) 

zeB{y,r) 

therefore, F(x,y, •) £ L^[B{y,r)^ due to the Holder inequality and our conditions on p. For 
z ^ B{y,r) due to the Holder inequality we have the estimate 

i/pi 

lK2Ur\mp, 



/ \F{x,y,z)\dz < ( 


[ \Ki{x,z)\^'dz] 


lx\B{y,r) \ 


. Jx\B{y,r) J 



and 

/ \Ki{x,z)\'^'dz < [ \Ki{x,z)\^'dz = I \Ki{x, z)\^' dz + [ \Ki{x, z)\^' dz , 

JX\B{y,r) Jx J B(x,r) Jx\B(x,r) 



where the first term on the right-hand side is finite due to the continuity of Ki , and the second 
one is estimated by (I2.1() . This proves the inclusion F{x,y, •) G L^{X). 

Now let Xo,yo G < r < i?, and x £ B{xo,r /2), y £ B(yQ,r/2), then 

\j{x,y) - J{xo,yo)\ < \F{x,y, z)\dz + \F{xo,yo, z)\dx 

■JB(yo,r) ■JB{yo,r) 

+ / \Fix,y,z)\dz + / \F{xo,yo,z)\dz 

Jx\B(v(^.R) Jx\B(va.R) 



IX\B{yo,R) JX\B(yo,R) 
JB(yQ,K)\B{yo,r) 



F{x, y, z) - F{xo, yo, z) dz. (2.8) 



Take e > and assume r < tq. For z G B{yQ,r) we estimate F(x,y,z) as in (|2.7|1 . then we get 
using Lemma 12 



/ \F{x,y,z)\dz < c sup Ki{x,y)\\W\\p / 

■JB{yo,r) xeB(xo,r), JB 



i{yo,r) x£B{xQ,r), 

y&B{yo,r) 

as r — > 0. On the other hand 



'B{yo,r) 



d{y, z)~ dz 



< Cr''-''-p = oil) 



/ \F{x,y,z)\dz < [Ki\p^^r[K2\p2,r\\Xx\B{xo,R)W\\p = o{l) as R ^ oo . 

Jx\B(xo,R) 

Finally, we conclude that r can be taken sufficiently small and R sufficiently large, such that the 
sum of the first four terms on the right-hand side of ()2.8() is less than e/2. Now it is sufficient 
to prove that at these fixed r and R the function 



/ 



F{x,y,z)dz 

R)\B{yo,r) 



is continuous as x G B(xo,r/2) and y G B{yQ,r/2). To do this, we note that with some C > 
the following estimate \F{x,y,z)\ < C'\W{z)\ takes place for all x G B{xo,r/2), y G B{yo,r/2), 
and z G B(yQ,R) \ B{yQ,r). Since W £ {B{yQ,R) \ B{yQ,r)^ , the requested continuity follows 
from the Lebesgue majorization theorem. □ 

As it was mentioned in the introduction, we are going to present the Green function in the 
form 

Ga,u{x, y; C) = Sa,u{x, y) + G'X^jix, y; (), 

where the second term must be continuous in X x X. Such a representation is trivial in the one- 
dimensional case: the Green function is continuous, and one can put Sa,u = 0. In dimensions 
v > 4 the problem makes no sense, as the following example shows: 

Example 7 (Four-dimensional Laplace operator). Consider the simplest case of the Lapla- 
cian in L^(M^). The Green function takes the form 

G{x, y- C) = . ^ . Ki (v^ \x - y\) , 
47r^|x — y\ 

where Ki is the modified Bessel function of the first order. Near the diagonal x = y one has 

^' = 4vr^|x-yP 8^ + ^' 
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with a continuous k. Therefore, for Ci)C2 £ res(— A), Ci 7^ (2, the difference 

G{x, y; Ci) - G{x, y; C2) ~ ^"^^^2^ log 1^ " ^1 

is a discontinuous function, so that the singularity cannot be chosen independent of the spectral 
parameter. 

Therefore, the only non-trivial cases remain u = 2 and = 3, which we will consider in the 
present article. 

Example 8 (On-diagonal singularity for the Laplace operator). Here we consider the 
case j4 = and [7 = 0, i.e. the case of the Laplace-Beltrami operator —A on the manifold 
X with u = 2 or V = 2>. Denote the Green function of —A by G{x,y;(^). Take y £ X and 
introduce polar coordinates {ry,uj), Vy = d{x,y), uj £ centered at y, then we have in a 

normal neighborhood Wy of y: 

\ Vy " OryJ OVy 

where the function 9y = 9y{ry,u;) is defined in such a way that in Wy, we have dx = 
ry~^6y{ry,uj)drydw. Since ry~^6{ry,uj) is the Jacobian for the inverse to the exponential map 

d 

in Wy, there holds 9y{<d,u>) > Cy > and —9y{0,uj) = for all lo G S'^"-^. Moreover, inf Cy > 
as y runs over a compact set in X. 
Denote now 

7;- log 



S{x,y) 



2vr d{x,yy 
1 



3, 



. 47r d{x,y) ' 

and for a fixed ( £ res(— A) denote K{x,y) := G{x,y;C,) — S{x,y). Then there holds 

rl9 f) 

(-A - C)K{.,y) = 9y'-^^S{;y) - CS{;y) =: L{x, y). (2.9) 

" Ory OVy 

It is clear that L{-,y) € LP'iWy), hence due to the Sobolev embedding theorem, x 1— > K{x,y) 
is continuous in Wy. Let us show that really K{x,y) is continuous in {x,y). To do this, we fix 
yo £ X and take tq > such that B(yQ, 2ro) C Wy^. We prove the following assertion: 

(CM) the map B{yQ,rQ) 9 y 1— > L{-,y) £ Lp'{B{yQ,rQ)) is continuous with respect to the norm 
topology of the space Lp'{B{yQ,rQ)). 

Let X £ C°°[X) such that suppx C B{yQ,2rQ), x{^) = 1 for x G B{yo,ro), and < x{^) ^ 1 
for all X £ X. Note that B{yQ,2rQ) is a normal neighborhood of y for all y £ B{yQ,2rQ), 
therefore we can assume that L{x,y) is defined for all x S X and y £ B{yQ,2rQ). Extend L 
by zero for y ^ B{yQ,2ro) and set T{x,y) = x{x)x{y)Li{x,y). It is clear that T £ K,cont{oi,p) 
where p is arbitrary number with 1 < p < 00, and a = 1 for = 3, a is any strictly positive 
number for u = 2. Using items (A) and (B) of Lemma El we can easily show that for every 

/ £ L'^{X) the mapping B{yQ,rQ) 3 y ^ / L{x,y)f{y) dy is continuous and the mapping 

B(,„,,,) 3 , ^ / ,1,0 „„,i„„„„3. This proves the assertion (CM). 

Returning to Eq. (|2.9() we see that K(-,y) tends to K{-,yo) with respect to the topology of 
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VFf (-B(yO)'"o))- Due to the Sobolev embedding theorem, this implies a uniform convergence in 
the ball B{yQ,r), i.e., 

lim sup \K{x,y) — K{x,yo)\ = 0. 

y^y^ xeB{yo,ro) 

This together with the continuity in x proves the required joint continuity in {x,y). Therefore, 
the functions S{x,y) are suitable on-diagonal singularities of the Laplace operator. 

Note that the proof of the separate continuity of the function K{x,y) is considerably simpler 
and can be found, e.g., in [21]. 

3 On-diagonal behavior for singular scalar potentials 

Below we will use the notation L^^^{X) = \Jg^pLl^^{X). 

Lemma 9 (Singularity is independent of the spectral parameter). Let u = 2 or 3, 

A G U G r{X), Ci,C2 G iMHa,u), then the difference y; Ci) - Ga,C7(x, y; C2) 

is continuous in X x X . 

Proof. The proof follows from the Hilbert resolvent identity for the kernels: Ra,u iCi) " 

RA,uiC2) = iCi - C2)RA,uiCi)RA,u{C2)- The integral kernel / GA,uix, z; Ci)GA,uiz,y; C2)dz 

Jx 

of RA,uiCi)RA,uiC2) is continuous due to LemmalHI^B). □ 

The previous lemma shows that for fixed A and U, the on-diagonal singularity in question 
exists; for example, as a singularity one can take GA,uix,y;Co) for a fixed Co G ^^s{Ha,u)- Our 
aim is to understand how the singularity depends on A and U. 

The following lemma shows that Green functions of Schrodinger operators with smooth 
potentials have the same on-diagonal singularity. 

Lemma 10 (Singularity for operator with smooth potentials). Let v = 2 or 3, A ^ 

{C°°{X)X , U,V £ V{X)nC°°{n), where n is a domain in X , then the difference GA,u{x,y;C)- 
GA,vix,y;C) has a continuous extension to all points {x,x), x £ In particular, ifVL = X, 
then GA,u{x,yX) - GA,v{x,y;0 S /Ccont(p) with arbitrary p>l. 

Proof. Fix a real E sufficiently close to —00 and take xq We show that in a neighborhood 
of {xq,xq) m. X X X, the difference y; E) = Ga,u{x, y, E) — Ga,v{x, y ; E) is the restriction 
of a continuous function in this neighborhood. Due to LemmalHIthe same will hold for all values 
of the spectral parameter. 

Let CIq be a bounded subdomain of $1 and contain xq; denote W = U + xnoiV it is clear 
that W E 'PiX). Since W — U is bounded with compact support, one has Ra,u{C) ~ RA,wiC) = 
Ra,u iOiW — U)RA,wiC)^ so that the difference 

GA,u{x,y;E)-GA,w{x,y;E)= [ Ga,u{x,z;E){W{z) -U{z))GA,w{z,y;E)dz 

Jx 

is continuous in X x X according to LemmaEfB). It remains to show that the function L{x, y) = 
GAy{x-, y; E) — Ga,w{x-, y, E) is continuous on JIq x f^o- To do this, let us note that in the sense 
of distributions the following equality holds: 

{{Ha,v). -E + Wl^y - E)L{x, y) 

= {W{x)-V{x))GA,w{x,y;E) + {W{y)-V{y))GA,w{x,y;E), (3.1) 
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where {Ha,v)x (respectively, {HA,v)y) means that Ha,v acts on the first (respectively, the 
second) argument in L; the bar means that we change the coefficients in H^y by the complex 
conjugate ones. The operator in the left-hand side of is elliptic in 0,q x 0,q with smooth 
coefficients, while the right-hand term vanishes in Oq ^ ^o- According to the elliptic regularity 
theorem L is continuous in f^o ^0- D 

The following Proposition contains our main result on the dependence of the on-diagonal 
singularity on singularities of the scalar potential. 

Proposition 11 (Preserving the on-diagonal singularity under singular perturba- 
tions). Let v = 2or?),A^ [C'^{X)Y , and Ui,U2 G 'PiX). If v = 3, assume additionally that 
U1 — U2& L'^^^{X). Then the difference GA,Ui{x,y; C) ~ G'a,i/2 (a^j y; C) ^-^ continuous in X x X 
for any ( S res{HA,Ui) n Tes{HA,U2)- 

Proof. For the sake of brevity we fix A and remove it from the notation, i.e. instead of Ga,u 
we will write Gjj etc. 

First of all, using Lemma ^ we choose functions Vi, V2 G C°°(X) semibounded below such 
that Wj := Uj — Vj = Yl^Li ^j,s, where Wj^g G L^^-" with 2 < pj^s < 00, s = 1, ... ,nj, j = 1, 2. 

For C G res(i7r/J n Tes{Hu2) the sets Vj := (Huj — C)C'^(X) are dense in L^(X), because 
Cq°{X) is an essential domain of both Hfj-^ and Hu^. As ip ^ T)j, one has 

Ru, (C)V' - Rv, (C)V' = Rv, {QW.Ru^ (C)^. (3.2) 

As the operators on the both sides of 1)3. 2|) are bounded and coincide on a dense subset, they co- 
incide everywhere, i.e. (|3.2j) holds for any G L?'{X). Combining Lemma|31and LemmaEl^B) we 
conclude that in the dimension two, the operator on the right-hand side of ()3.2|) has a continuous 
integral kernel, which together with Lemma [Till implies the conclusion of the proposition. 

Let us consider the dimension three more carefully. To be shorter, we remove the dependence 
of the resolvents on C, from the notation. We have the following chain of equalities: 

Ru^ - = Rvi - Rv2 + Rvx^iRui - RV2W2RU2 

= Rv, - Rv2 + Rv.WiRu, + Rv2W2{Ru^ - RU2) - RV2W2RU, 

= Rv, - Rv2 + Rv2W2{Ru, - RU2) + Rv2{Wi - W2)Ru, + {Rv, - Rv2)WiRu,. 

Therefore, (1 - Rv2W2){Ru^ - RU2) =:L = A + B + C, where A := Ry^ - Rv2, B := Rv2iWi - 
W2)Ru^, C := {Rv, - Rv2)WiRu,. 

Due to Lemma [TUl the operator A has an integral kernel from /Ccont(p) with arbitrary p, 
p > 1. Since Wi — W2 S L\^^{X), the operator B has an integral kernel from /Ccont(oo) due to 
Theorem 1^ and the items (A), (B) of Lemma IHl As Rv2 — Rvi £ ^cont(p) with arbitrary p> 1 
(Lemma llfl|). the integral kernel for C is from /Ccont(oo) due to Theorem |2l again and the items 
(A), (C) of Lemma ini Therefore, the operator L has an integral kernel L{x,y) = L{x,y;() G 
^cont(c>o). Now we note that the multiplication by W2,s is a continuous mapping from L°°{X) 
to LP2,s(|x). At the same time, as £ ^^00111(15^)5 P > 1, the resolvent Rv2 is a bounded 
operator from each L^'2.«(X) to L°°(X) due to Lemma IHl Since L = (1 — Rv2W2){Rui — Ru2)i 
we can combine Theorem |21 and Lemma El to show that the operator L is a bounded map from 
LP{X) to L°°{X) for any p with 3/2 < p < oo. Since \L(x, y;C)\ = \L{y, x;C)\, we see from (|2.3p 
that L{x,y) G /Ccont('7) for any q with 1 < g < 3. 

One can find C such that ||-Rv2 (C)^2||oo,oo =■ a < 1 (see [20]), therefore, the operator 
1 — RV2W2 acting in L°°{X) is invertible and for any n G N there holds 

n-l 

Ru, - Ru2 = Y,(^V2W2)''L + (1 - Rv2W2)-\Rv2W2rL. (3.3) 

k=0 
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Applying iteratively Lemmas El and IHIJ A) and taking into account Theorem |21 we can show that 
the operators {Rv2W2)'^Rv2 have integral kernels from /C cont(/?fe5 oo) with < 1. At the same 
time, all these operators are bounded from to L°°{X) for any p with 3/2 < p < oo. Using 

the same arguments as for L above, we conclude that these kernels are in fC contl/^fc) q) for any q 
with 1 < g < 3. Applying now Lemma (C) one proves that the first term on the right-hand 
side has a continuous integral kernel. 

Denote T„ := {1 — Rv2W2)~^ {Rv2W2)^~^ Rv2y this operator is bounded from each LP^-'^X) to 
due to the Gelfand-Dunford- Pettis theorem, this is an integral operator with an integral 
kernel Tn{x,y). The second term in (|3.3|) takes the form TnW2L, and by virtue of Lemma |31 

this is also an integral operator with the kernel Sn{x,y) := / Tn{x, z)W2{z)L{z,y) dz. Prom 

Jx 

the other side, one can write Sn{x,y) = TnW2ly{x), where ly{x) := L{x,y). Note that for each 
y € X there holds ly € L°°{X), and the operator TnW2 is a bounded mapping from L°°(X) to 
L^{X) with the norm ||r„VF2||oo,oo < ||(1 - «V2W^2)-'|L,^ • Py^VFsir^,^ < a^/{l - a). 

Now let us fix xo S and take a bounded open neighborhood 17 of xq- It is clear that 
ll^ylloo < Co for all y € $7 with a certain cq > 0. Therefore sup^ y^Q \Sn{x,y; ()\ < CQa"'/(l — q). 
Take e > and choose n such that cna"/(l — a) < e. Prom Eq. ()3.3|1 we have in J7 x the 
relation Gui{x,yX) — Gu2{x,y',C) = Kn{x,y) + Tn{x,y), where is continuous and \Sn\ < £• 
As e is arbitrary, this means that Gfj^ (x, y; () — Gu2 {x^ y; C) is continuous in 17 x J]. Since xq £ X 
is arbitrary, the lemma is proven. Due to LemmalHl this holds for all C £ Tes{Hy^)nres{Hv2) ■ D 

The following example shows that the condition Ui — U2 € L^^{X) can not be omitted in 
dimension three. 

Example 12 (Coulomb potential in three dimensions). Let X = M.^, A = 0, and U = 
q/\x\, i.e. H = Ha,u = —/S.+q/\x\. Clearly, U ^ Lj|^^(R'^). The Green function can be calculated 
explicitly [22]: 



r(i 



47r|x — y\ 



where := |x| + \y\ + \x — y\, r] := \x\ + \y\ — \x — y\, k = —q/\^—AC, -Wk,i/2 and W«;_i/2 are the 
Whittaker functions, 

M,,i/2(x) = e^'/2x$(a, 2; x) , W^^y2{x) = e'/^x^{a, 2; x) . (3.5) 

Here $(a, c;x) and ^{a,c;x) are the Kummer function and the Tricomi function, respectively. 
We prove in Appendix^ the asymptotics 

G(x,0;C) = j^ + ^log|x|-^ 
4-k\x\ Att 4tt 

+ £ (V' (l + + log + log(2/e) + 2Ce^ + 0{\x\ log |x|). (3.6) 

Therefore, the singularity for G{x, y; Q contains an unavoidable logarithmic term and is different 
from the standard three-dimensional singularity. 



4 Dependence of the singularity on the magnetic field 

Lemma 13 (Singularity due to the magnetic field in two dimensions). Let v = 2, 

then for any A £ [C°^{X)Y the difference GA,oix, y; C) — ^0,0(2;, y, C) 'is continuous in X x X if 
C G res{HAfl) n res(i?o,o)- 
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Proof. Let xq be an arbitrary point of X. We show that the difference GA,oix, y, C)—Gofi{x, y; C) 
is continuous in a neighborhood of {xq,xo) for at least one value of the spectral parameter 
due to Lemma IHl this difference is continuous for all admissible spectral parameters. 

Take two sufficiently small numbers r and tq with < r < tq. Fix a function (p £ Cq°{X) 
such that supp(/> C B{xQ,ro), (j){x) = 1 as x E B{xQ,r). Denote for brevity Hq := -ffo,o> 
Hi := Ha,o, H2 := H^Afl] the corresponding Green functions will be denoted by Go, Gi, and 
G2, respectively. 

In B{xQ,r) X B(xq, r) for real C, sufficiently close to —00 one has in the sense of distributions 



(((i^i), - C) + {{H2)y - C)) (Gi(x, y- - G2(x, y- C)) 



0, 



therefore, due to the elliptic regularity, the difference Gi{x,yX) — G2{x,y;C) is continuous in 
B{xo, r) X B(xo,r). Now we are going to show that G2{x, y; C) — Go{x, y; C) is continuous. Since 
Hq and H2 are uniformly elliptic operators with G'^-bounded coefficients, we are able to use 
estimates for the Green functions and their derivatives obtained in [19]. First of all, 



Go(x,y;C), G2(x,y;C) G /Ccont(A, 



(4.1) 



for arbitrary A > and q G [1,cxd] (see Theorem I^J. Moreover, for C close to —00 both these 
kernels are smooth outside the diagonal x = y, and according to [19, Theorem Al.3.7] we have 



dxGoix,y;C) 



< G 1 + 



log d{x,y)\ 
d{x,y) 



-wd{x,y) 



J 



1,2, 



where d is any first order derivative taken in canonical coordinates, and G,uj > 0. Additionally, 
by [19, Theorem Al.2.3] for any p > 1 there exist e, G' > such that 



sup 



dxGo{x,y;C) 



I d(x,y)>r 

This implies the inclusion 



p 

e'"^('^^y)dy + sup 

y Jd{x,y)>r 



dxGo{x,y;C) 



dxGo{x,y;C) S ^cont(l + A,( 



e''ii^'y)dx<C', J = 1,2. 



(4.2) 



with the same A and q as in 1)4. 

In canonical coordinates in B(xQ,rQ) both Hq and H2 are given by symmetric second-order 
elliptic expressions with the same principal symbol, in particular, the difference T := H2 — Hq 
is defined by a first order differential expression, T = 61 (x) di + 62(2;) 82 + c(x), where 61, 61, 
c are compactly supported smooth functions. For the functions of the form ip = {Hq — C)4' 
with (f) G C^{X) we have {H2 - C)0 = {Hq + T - C)i?o(C)V' = (l + TRq{C))'iP, therefore, 
Ro{C)'4' ~ -^2(0^ = R2{C)T Ro{C)ip ■ In terms of integral kernels this means 



GQ{x,y;C)i;{y)dy - / G2{x,y; C)ip{y)dy 



X 



X 



X 

G2{x,z;C) 
X Jx 



G2{x,z;C) 



bi{z) di + b2{z) 82 + c{z) 



X 



Go{z,yX)4'{y) dy dz 



bi{z)Ki{z,y;C) + b2{z)K2{z,y-C) + c{z)Go{z,y;C) 



iP{y)dydz, (4.3) 



where 

Ki{z,y;C) ■■= dziGo{z,yX), K2{z,yX) := dz2GQ{z,y, ()■ 

According to the general theory of elliptic operators, the set {Hq — ()C^{X) is dense in all 
L'P{X) with any p, 1 < p < 00, if C is sufficiently close to —00 [19, Section Al.2]. Due to the 



11 



estimates (|4.1() . (|4.'2() . and Lemma [SJ the kernels Ki and define bounded operators from 
L'^{X) to L°°(X) for arbitrary (7 > 2; denote these operators by Ki{C) and i^2(C)- this 
notation, the expression in the right-hand side of (|4.3|) can be rewritten as 

i?o(C)V' - ^2(0^- = [R2{C)biKi{C) + R2{C)b2K2{C) + ^2(0 cRiC)]ip. 

The operators in the both sides are bounded from L'^{X) to L°°{X) with any q > 2 and coincide 
on a dense subset, therefore, the corresponding kernels coincide, i. e. 

Go{x,y;C) - G2{x,y;C) = / G2{x,z;C)bi{z)Ki{z,y;C)dz 



X 

+ f G2{x,z;C)b2{z)K2{z,y;C)dz+ [ G2{x,z;C)c{z)Go{z,y;C)dz. (4.4) 
Jx Jx 

By Lemma ini(B), the function on the right-hand side of (|4.4|) is continuous. □ 

The three-dimensional analogue of Lemma [THl is not true as the following example shows. 

Example 14 (Three-dimensional Landau Hamiltonian). Consider in L^(M'^) the vector 
potential of a non-zero uniform magnetic field. By a suitable choice of coordinates one can 
assume that the field is directed along the xs-axis, i.e. the magnetic strength vector is B = 
(0, 0, 27r^j;3), where ^ > is the density of the magnetic flux through the plane {xi,X2)- Choose 
the symmetric gauge for the the magnetic vector potential, A(x) = x x, then H := Ha,o 
takes the form 

rr f- 9 . \2 /. 5 ^ \2 8^ 

H = U- <X2 + hTT^Xi) -T-2' 

V oxi J V 0x2 / dx\ 

and the corresponding Green function is G(x,y; Q = <I>(x,y)F(x — y; C), where 

-■oo exp[-7r|^|(x^(e* - 1)"^ +x^t-ii 



i^(x;C) 



(1 — e *) exp 



1 c 



dt , (4.5) 



.2 47r|^|, 

x_L = (xi,X2,0) and xy = (0,0, X3) [23]. In Appendix IHl we prove the asymptotics 

G(x, y; C) = -T-r--T + 7 ^ Z o ; o - 7dl7 + - y|) (4.6) 



47r|x - y| 4 V vr y \2' 2 47r|^| 

as |x — y| — > 0; here Z{z;u) is the generalized Riemann (^-function (also known as the Hurwitz 
function) . Therefore, the on-diagonal asymptotics is 

^(^^y) = -n r = t~\ r 

47r|x — y| 47r|x — y\ 



5 Summary of results 

We summarize some corollaries from the proven assertions in the following two theorems. 

Theorem 15 (On-diagonal singularities of the Green functions in dimension two). 

On a two-dimensional manifold of bounded geometry X, for any vector potential A G [C°°(A)]2 
and scalar potential U £ V{X), the Green function Ga.u of the Schrddinger operator Ha,u = 
— -|- U has the same on-diagonal singularity as that for the Laplace- Beltrami operator, i.e. 

Ga,u{x, C) = ^ log + G'Xjjix, y; (), 

where G^^^ is continuous on X x X . 
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Proof. Proposition 111! shows that the singularity does not depend on the scalar potential U G 
V{X), and Lemma [T^ shows that it is independent of the magnetic potential. Therefore, the 
singularity coincides with that for the Laplacian, see Example |S1 □ 

Theorem 16 (On-diagonal singularities of the Green functions in dimension three). 

Let X be a three-dimensional manifold of bounded geometry. For U G ^{X) and A G [C°^{X)]^ 
consider the Schrddinger operator Ha,u = —^A + U and its Green function GA,u{x,y]C)- If 
Ui,U2 G 'P{X) and Ui — U2 € L'I^^{X), then the Green functions Ga,Ui Ga,U2 have the 
same on-diagonal singularity {i.e. Ga,Ui ~ Ga,U2 ^-^ continuous in X x X). In particular, for 

loc( 



any U G n Lf;^^{X) there holds 



Go,u{x,y;C) = ^^^[^ + G^^^(x,y;C), (5.1) 
where G^Q^ is continuous in X x X . 

Proof. The theorem is a simple corollary of Proposition II 11 and the formula 1)5. 1(1 follows from 
Example El □ 

Remark 17. Contrary to the two-dimensional case, the singular term of the Green function for 
the three-dimensional Schrodinger operator Ha,u does depend on the scalar potential U as well 
as on the magnetic vector potential A. In particular, if A is the vector potential of a uniform 
magnetic field B in X = R^, then instead of ()5.1|) we have 

GA,o(x,y;C) = r exp (!^(^^) + G'Xfii^,y;C) , 

4:Tt\x — y\ \ 1 J 

see Example II 41 On the other hand, the dependence on scalar potentials is shown in Example ll2l 
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A OfF-diagonal asymptotics for the inverse harmonic 
oscillator in dimension one 

The Green function G{x,y\C,) for the inverse harmonic oscillator H = —dP/dx^ — uP'x^ j^.., has 
the form 

G(x,y;C) = 

xU{- iC/iv, e-'^l^u^l'^ max(a;, y)) xU{- iC/iv, e'^^^^J"!'^ max(-x, -y)) , (A.l) 

where 9C > ^^id U{a,x) is the Weber function, see [24, Chapter 19]. Using [24, no. 19.8.1], 
for large z one obtans U{a,z) = e~^^^^z~2~"' u{z), where lim^^oo ^^(-z) = 1- Returning to the 
Green function we see that for fixed x and large y one has (assuming y > x) 

Gix, y; = U{ - ^C/u;, -e^/^cV^ ^) ^ 

V2TTUJ (e-W4t^l/2y)-</'^+2 

where l[my-,oov{y) 7^ 0. Therefore, for large \x — y\ the Green function has only a polynomial 
decaying. 
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B On-diagonal singularity for the Coulomb Hamilto- 
nian 

Here we prove the asymptotics ()c{.6|) . 

We are interested in asymptotics of the functions x G{x,xq;C) as x — > xq at fixed 
C G res(ff) and xq £ M^. As the potential is smooth outside the origin, the Green function 
has the standard on-diagonal asymptotics if xq ^ 0. We consider the case xq = 0. We have 
A^K,i/2(0) = 0, M^_^/2(0) = 1' therefore, 

G{x,0;O = ^^^^VF.,i/2(2v^|x|). 

Consider the following expansions (cf. items 6.1(1) and 6.8(13) in [25]): 

a a(a + 1) 9 
^a,2;x) = 1 + -X + ^-^x^ + . . . , 



^) = ^fR + ^) ^ + g ^ r(a)(. + l)!fe! 

A^ix~^ + ^0 + Aix + ^22;^ + . . . + -Bo log X + Bix log x + B2X^ log x + 



where 



1 _ ^(a)-^(i)-^(2) _ fl(^(fl + l)-V,(2)-V;(3)) 

^-'"i>)' I>^T^ ' 2r(a-l) ^ 



a(a + l)(V>(a + 2)-V;(3)-V>(4)) 1 „ « „ a(a + 1) 

^2 = TTTw T\ ' -°0 = — -Dfl = -r-^, 02 — 



12r(a-l) ' " r(a-l)' " 2r(a-l)' ' 12r(a-l) 

Using (|3.5|) . we get 



W^^il2{x) =^-1+ (^Ao - X + So2;logx + O(|x2logx| 



1 /V'(a) -V'(l) - ^^(2) 1 \ 1 , ^,,2, 

+ ^ ^ ^ / y ^ _^ _^ -xlogx + 0(|x^logx|) 



r(a) V r(a-i) 2r(a)y r(a-i)' 

Since -0(1) = —Ce-, "0(2) = 1 — C'^;, where Ce is the Euler constant, we get 1)3. 6(1 after some 
trivial algebra. 



C On-diagonal singularity of the three-dimensional 
Landau Hamiltonian 

In this appendix, we are going to prove the asymptotics (|4.6() . 

Set in the integral (|4.5|) x_l = and denote xy = z. Then after the change of variables t ^ 
in this integral, we obtain 

G(0, 0, z; 0, 0, 0; C) = ^ / ^ dt , (C.l) 

1 C 

where a = 7r|^| and c = — — -. Represent now G(0, 0, z; 0, 0, 0; C) = fi{z] Q + f2{z; (), where 

2 47r|^| 



2tt Jq Vl-e-* 
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Changing the variable t ^ t ^ and using the relation 

/ expi-bt"^ - c/t^) dt = -(vr/6)i/2 expi-2{bc)^/^) 
Jo 2 

(see [26], V. I, Formula 2.3.16.3), we obtain fi{z;C) = exp ( - (27r|^| - C)^/^kl)/(47r|z|) , or 
G(0,0,2;0,0,0;C) = (47r|z|)~^ + ^(z; C), where 

5(^;C) = -^(2vr|e|-C)'/' + /2(^;C). (C.3) 

It is clear that the function g is continuous with respect to z and analytic with respect to 
C S ices{HA,o)- We can rewrite HC.1|) in the form 

^ r exp(-vr|g|z^t-^) 1 

(l-e-t)exp^i-4^)tjVt 



Let = (e* — 1) ""^ — t ^; the function h is real-analytic on the whole line, h(t) ^ as t ^ +cxd 
and h(t) — > —1 as t ^ — oo. Therefore, h is bounded on R. Let us represent F{x; (") in the form 



+ 



F{x;C)-- 
exp(-7r|^|x2t"i 



exp(-7r|^|x^t-i) 



(1 — e *) exp 



1 



c 



■dt 



Vt 



(1 — e *) exp 



c 



1 



2 47r|^| 



Vt 



2 47r|^|, 

{exp[-7r|e|xx/i(t)] - 1} = /i(x, C) + /2(x, C)- 



(C.5) 

It is easy to show that I2 is a continuous function in the domain x G M'^, Re^ < 27r|,^|. Let 
us show that /2(x, C) locally uniformly with respect to ReC < 27r|^|, as x — > 0. It is 
sufficient to show that 



^(x,C) 



exp(-7r|^|x2t-i) 



(1 — e *) exp 



c 



1 



2 47r|^| 



Vt 



|exp[-7r|^|xi/i(t)] -l\ dt^O 



locally uniformly with respect to C G M, C < 27r|^|, as x — > 0. Fix C £ M, < 27r|^|. Since 
x^ < x^, we have |exp[— 7r|^|x^/i(t)] — l| < constx^ in a neighborhood of the point (0,0, z). 
Therefore, using (|C.4() . we get 



^(x, C) < cx^ 



exp(-7r|^|x^t-i) 



cx 



° (l-e-*)exp 



1 



c 



Vt 



,2 47r|e|, 

and we get the required limit. Using ()C.4|) again, we obtain 

1 



^i(|x|,C) 



|e|V2|x| 



+ /(|x|,C)- 



(C.6) 



Using (|C.5|) and (|C.6() we get 



n/ ^\ 1 exp [7r«^(x_L A y_L)] , ^, 
G(x, y ; C) = , , + F(x, y; C) , 



4tt 



|x- y| 
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where F(x, y;^) is jointly continuous with respect to (x, y) G x for aU ^ G res(ffA,o)- 

Denote Q{C) = hm|x-y|_»o -^(X) Yi C) j this hmit is independent of x and y since -F(x, y; is 
invariant with respect to magnetic translations T^, a G M'^: ra/(x) = exp[7ri^(a_|_ Ax_|_)] /(x — a). 
From (|4.5() we obtain 



Aq(0 = ^— 



exp 



c 



4vr|e| 



(1 



/■oo 

Using Equation (1.10.4) from [25] we get / t 

Jo 

obvious relation dZ(s,v)/dv = —sZ{s + l,v) implies immediately 



r(s) Z(s, v) and the 



QiO 



1/2 



1 1 



c 



2 ' 2 47r|^| 



with a constant C G M. To determine C we compare HC.7|) with ()C.3|) in the limit 
Since Q(C) = 5(0; C), we have from (ITThIi and (ITT^ : 



(C.7) 



-oo. 



g(C)-^(2vr|e|-C)'/' 



as ~^ "OO . 



On the other hand, by the Hermite relation (see (1.10.7) from [25]) there holds Z(l/2,r;) + 
2tii/2 _> as JRv +00. Comparing the two last relations with (|C.7() . we get C = 0. Thus, 
()4.6|) is proven. Note that the expression for Q(C) was obtained at the physical level of rigor 
in [27] and can be found also in [7]. 
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